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We unearth the interconnection between various analytical methods which are 
widely used in the current literature to identify integrable nonlinear dynamical sys¬ 
tems described by third-order nonlinear ordinary differentiable equations (ODEs). 
We establish an important interconnection between extended Prelle-Singer proce¬ 
dure and A-symmetries approach applicable to third-order ODEs to bring out the 
various linkages associated with these different techniques. By establishing this in¬ 
terconnection we demonstrate that given any one of the quantities as a starting 
point in the family consisting of Jacobi last multipliers, Darboux polynomials, Lie 
point symmetries, adjoint-symmetries, A-symmetries, integrating factors and null 
forms one can derive the rest of the quantities in this family in a straightforward and 
unambiguous manner. We also illustrate our findings with three specific examples. 
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1. Introduction 

In a previous paper (Mohanasubha et al. Proc. R. Soc. A 2014) we have established 
a connection between extended Prelle-Singer procedure with five other analytical 
methods which are widely used to identify integrable systems described by second- 
order ODEs and brought out the interconnections between Lie point symmetries, A- 
symmetries, Darboux polynomials, Prelle-Singer procedure, adjoint-symmetries and 
Jacobi last multiplier. We have also illustrated the interconnections by considering 
a nonlinear oscillator equation as an example. A natural question which arose was 
whether these interconnections exist in higher order ODEs as well. In this paper 
we consider this problem for third-order ODEs and come out with certain new and 
interesting results. 

We start our investigation with the following question. Given any one of the 
quantities as a starting point in the family, consisting of multipliers, Darboux poly¬ 
nomials, Lie point symmetries, adjoint-symmetries, A-symmetries, integrating fac¬ 
tors and null forms can one obtain the rest of all quantities in this family? In the 
case of third-order ODEs exploring the interconnections between the above said 
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quantities is difficult both conceptually as well as technically which will be clear 
as we proceed. In fact, we get an affirmative answer to our questions. We divide 
our analysis into two parts. In the first part, as we did in the case of second-order 
ODEs we try to interlink the quantities by introducing suitable transformations in 
the extended Prelle-Singer procedure. Once the interconnections are identified, in 
the second part, we demonstrate that starting from any one of the quantities in this 
family one can derive all the other quantities in a simple and unambiguous way. 

In the first part, the interplay between various methods is brought out by es¬ 
tablishing a road map between extended the Prelle-Singer procedure and the other 
methods, namely (i) Jacobi last multipliers, (ii) Darboux polynomials, (iii) Lie point 
symmetries, (iv) adjoint-symmetries and (v) A-symmetries. In the above we have 
listed out the methods in chronological order. In other words JLM was the one 
invented first and the last one in this direction is the extended Prelle-Singer pro¬ 
cedure. So naturally we try to accommodate other methods into extended Prelle- 
Singer procedure. In the later procedure one has three equations to integrate to 
obtain the null forms U and S and the integrating factor R. To interconnect these 
quantities we introduce three new variables V, X and E, which are to be determined 
through the transformations, U = S = y and i? = ^ so that one can 

rewrite the determining equations for U, S and R in terms of X, V and E, respec¬ 
tively. Doing so, we find that the determining equation for the function F exactly 
coincides with the determining equation for the Darboux polynomials which in turn 
establishes a connection between the integrating factors and Darboux polynomials. 
Now properly combining the expressions which arise in the extended Prelle-Singer 
procedure with the A-symmetry determining equation we obtain the following re¬ 
lation between U and S with A (see Sec.3 for details), namely U = _ _ 

This interlink is vital and new to the literature. Since the null forms U and S are 
now expressed in terms of V and X, this interlink can also be written in terms 
of V and A. Doing so we obtain an equation solely in terms of V and A, that is 
D^\y\ — D\V]{(j)x -I- A) -b IA(iA[A] -b A^ — ())*) = 0 which in turn determines A if E is 
known or vice versa (Here x, x, x) is the function defining the third-order ODE, 
see Eg. (12.II) below). We assume that A is known. In this case V can be determined 
by solving this equation and from the latter X can be fixed and so U, S and R 
can be determined. The Darboux polynomials can be determined either from the 
multipliers or from their own determining equations. This in turn establishes the 
interconnection between all these methods. One important result which we observe 
in this part is that for a given A, the null forms U and S are not unique which can 
be seen from the above relation. Because of this, one may find multiple null forms 
for a given equation and one has to choose the correct form of U and S in order 
to obtain an independent integral. From the remaining forms of U and S one may 
obtain other useful information such as symmetries, integrating factors, Darboux 
polynomials and Jacobi last multiplier; however the resulting integrals turn out to 
be dependent ones. This feature is different from second-order ODEs where one has 
one-to-one correspondence between the null forms and A-symmetries (A = —S). We 
illustrate this important point by considering three examples. 

Next we move on to the question that given any one of the system quantities how 
to derive the remaining quantities in this family. Using the interconnections which 
we identified above, and starting from any one of the quantities, we establish road 
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maps to connect to all the other quantities. After a detailed analysis we find that it 
is sufficient to consider any one of the following three cases as starters, namely (i) 
Lie point symmetries, (ii) integrating factors or (iii) Darboux polynomials. Then 
the remaining quantities can be uniquely determined. The other three possibilities, 
namely (iv) A-symmetries, (v) null forms and (vi) JLM are only subcases of the 
previous three cases, respectively. The exact road map in all the cases is explicitly 
demonstrated in Sec.4. To illustrate the ideas introduced in this paper we consider 
three examples. 

Since some of the methods applicable for the third-order ODEs are not much 
discussed in the literature as in the case of second-order ODEs, to begin with we 
review briefly the six different analytical methods in the next section (Sec.2) in 
order to be self contained. 


2. Methods 

In this section we briefly recall the six analytical methods with reference to the 
third-order ODEs which are widely used in the contemporary literature to derive 
one or more of the following quantities, namely (i) multipliers, (ii) integrating fac¬ 
tors, (iii) symmetries and (iv) integrals. We present the methods in chronological 
order. 


(a) Jacobi last multiplier method 

The Jacobi last multiplier method was introduced by Jacobi in the year 1844 (Ja¬ 
cobi, 1844,1886), but laid dormant for a long time. Nucci et al. have demonstrated 
the applicability of this method in exploring non-standard Lagrangians associated 
with certain second-order nonlinear ODEs (Nucci & Leach, 2008). The attractive 
feature of the method is that if we know two independent Jacobi last multipliers 
Ml and M 2 of the given equation, then their ratio yields a first integral for the 
given equation. 

We consider a third-order ODE, 


'x = (j)(t, X, X, x), 


( 2 . 1 ) 


where (/) is a function of t, x, x and x. We rewrite the third-order ODE (EH) equiv¬ 
alently as a system of following three first-order ODEs in the absence of explicit 
t-dependence 

i = f{x,y,z), y = g{x,y,z), z = h{x,y,z), (2.2) 

where /, g and h are suitably chosen functions of x, y and z. The analysis can 
be extended in principle to the general case with t-dependence included, though 
it becomes more involved now (see for example, Clebsch, 2009). We consider two 
integrals Ii{x,y,z) and l 2 {x,y,z) of (|2.2I) whose total differentials are given by 


.dh .dh .dh „ .dh .dh .dh „ 

X— - Vy— + z— = 0, X— - \-y—+z— = 0. 

ox oy oz OX Oy Oz 


(2.3) 
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The equations of motion which satisfy these two conditions can be written as 

= fi{x,y,z), 


X = 


y = 


dll 9/2 dll 9/2 

dy dz dz dy 

_ 9/1 dh 

dz dx dx dz 

\ dx dy dy dx 


= h{x,y,z), 
= f 3 {x,y,z), 


(2.4) 


where we have defined for simplicity the expressions inside the parenthesis as 
fi(x,y,z), f 2 {x,y,z) and f 3 {x,y,z). We assume that these equations be the ones 
coming out from the original equation (IQ) after multiplying by an integrating 
factor M{x,y,z). Then comparing Eas. (j2.2l) and (12.4L we have 


fi = Mf{x,y,z), f 2 = Mg{x,y,z), fs = Mh{x,y, z). 
Interestingly one may also observe from (EH) that 

dfi , 9/2 , 9/3 


dx 


dy dz 


(2.5) 


( 2 . 6 ) 


Substituting (j2.5p in (1^ and expanding it we can get the following equation for 
the multiplier M, namely 


9M 


dM 


dM 


fix, y,z)—+ gix, y,z)— + hix, y, z)—+[— + — + —] M = 0. (2.7) 


dx ' ^ dy 

Ea. (l2.7|) can be simplified to yield 


dz 


df dg dh 


dx dy dz 


where the differential operator D = x^ + y-^ + i^. Choosing fix,y,z) = y, 
gix, y,z) = z and hix, y, z) = (p, Ea. (l2.8|) can be written as 

D[\ogM] + p^ = 0. (2.9) 

Substituting the given equation in (EH) and solving the resultant equation one can 
obtain the multiplier associated with a given third-order ODE (Jacobi, 1844,1886; 
Nucci & Leach, 2008). Note that there can be more than one multiplier for a given 
ODE ()2.1I) . We further remark here that when the function (p in (12.11) also de¬ 
pends on t explicitly, the above operator D will be replaced by the total differential 
operator D = ^-^ + x^ + y-§^ + z-§^.^ 

An important application of multipliers is that one can determine the integrals 
associated with the given equation by evaluating their ratios. If Mi and M 2 are two 
multipliers, then it is straightforward to check from (1^ that one can identify an 
integral as I = ^ (Nucci, 2005). So if we have sufficient number of multipliers we 
can obtain the necessary integrals to prove the integrability of EH- 
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( 6 ) Darboux polynomials approach 

Darboux polynomials method was developed by Darboux in the year 1878 (Dar¬ 
boux, 1878). It provides a strategy to find first integrals. Darboux showed that if 
we have -I- 2 Darboux polynomials, where n is the order of the given equa¬ 

tion, then there exists a rational first integral which can be expressed in terms 
of these polynomials. In the case of third-order DDEs, we can get eight Darboux 
polynomials for a given equation (Darboux, 1878). 

The Darboux polynomials determining equation is given by the following ex¬ 
pression 

D[F]= g{t,x,x,x)F, (2.10) 

where D is the total differential operator and g{t, x, x, x) is the cofactor (Darboux, 
1878). Note also that for a given F, and an integral I of (12.ip . the quantity f{I)F, 
where / is arbitrary, is also a solution of (I2.10p for the same cofactor. Using this 
fact and solving Ea. (l2.10p . we can obtain the Darboux polynomials (F) and the 
cofactors (g), see Mohanasubha et al. (2014a) for further details on the method. 


(c) Lie symmetry analysis 

Lie symmetry analysis is one of the powerful methods to investigate the inte- 
grability property of the given ODE of any order n where one first explores the 
symmetry vector fields associated with it. The Lie symmetry vector fields can then 
be used to derive integrating factors, conserved quantities and so on (Olver 1993). 

We consider a third-order ODE of the form (EH). The invariance of Eq. E3D un¬ 
der an one parameter group of Lie point symmetries, corresponding to infinitesimal 
transformations. 


T = t + e^{t,x), X = X + er]{t,x), e <C 1, (2-11) 


where ^{t,x) and ri{t,x) are the coefficient functions of the associated generator of 
an infinitesimal transformation and e is a small parameter, is specified by 


^ dt 


dfP , (l)i9(/> (2)90 (3) „ 


( 2 . 12 ) 


Here r/9) ^ r)( 2 ) and are the first, second and third prolongations, respectively, 
of the infinitesimal point transformations ( 12 . 111 ) and are defined to be 

7^(1) = fi — x^, — x^, — x'i. (2.13) 


In the above over dot denotes total differentiation with respect to t. Substituting 
the known expression 0 in ( 12 . 121 ) and solving the resultant equation we can get the 
Lie point symmetries associated with the given third-order ODE. The maximum 
number of admissible Lie point symmetries for a third-order ODE (12.11) is seven 
(Olver 1993). The associated vector field is given by D = ^^ -|- 77 ^. 

One may also introduce a characteristics Q = rj — x^ and rewrite the invariance 
condition (12.1211 in terms of a single variable Q in the form 

D3[Q]=0iD"[g] + 0iD[Q]+0,Q. (2.14) 

One can deduce the coefficient functions ^ and g associated with the Lie point 
symmetries from out of the class of solutions to (I2.14p which depends only on x 
and t, and also has a linear dependence in x. 
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(d) Adjoint-symmetries 

In general, for systems of one or more ODEs, an integrating factor is a set of 
functions, multiplying each of the ODEs, which yields a first integral. If the system 
is self-adjoint (that is the adjoint of the linearized symmetry condition is the same 
as the linearized symmetry condition) then its integrating factors are necessarily 
solutions of its linearized system. Such solutions are also the symmetries of the given 
system of ODEs. If a given system of ODEs is not self-adjoint, then its integrating 
factors are necessarily solutions of the adjoint system of its linearized system. These 
solutions are known as adjoint-symmetries of the given system of ODEs. 

The adjoint ODE of linearized symmetry condition (12.141) can be written as 
(Bluman & Anco, 2002) 

L*[x]K{t,x,x,x) = D^[K\+D‘^[(j)^N\ - D[(j)iK\ + (j.K = 0. (2.15) 

where D is the total derivative operator which is given by D = 

Solutions of the above equation are known as adjoint-symmetries. The adjoint- 
symmetry of the Eq. (EH) becomes an integrating factor of (EH if and only if 
A(t, X, x) satisfies the adjoint invariance condition 

Ata + ^xxi + 2Aj, -|- xA±x + {4>^)xx = 0. (2.16) 

Once we know the integrating factors, we can find the first integrals. Multiplying 
the given equation by these integrating factors, and rewriting the resultant equation 
as a perfect differentiable function, 

A{t,x,x,x){x - (l}{t,x,x,x)) = ^{I), (2-17) 

we can identify the first integrals. 

(e) symmetries 

All the nonlinear ODEs do not necessarily admit Lie point symmetries. Under 
such a circumstance one may look for generalized symmetries (other than Lie point 
symmetries) associated with the given equation. One such generalized symmetry is 
the A-symmetry. Let V = ^{t,x)^ ri{t,x)-^ be a A-symmetry of the given ODE 

for some function A = X{t,x,x,x). The invariance of the given ODE (12.11) under 
A-symmetry vector field is given by 1/D>(3)]^'^' _ cj){t,x,x,x)) — 0, where 1/D>(3)] 
is given by -I- r]-§^ -§^ + Here 771 '^di)]^ 7^D.(2)] 

and are first, second and third A- prolongations respectively whose explicit 

expressions are given by (Muriel & Romero, 2001, 2008, 2009) 

y^[A.(i)] ^ {D + \)q{t,x) - {D + X){£_{t,x))x, 

^[A.(2)] ^ (D + A)?7[^’(i)l(t,x,i)-(D + A)(^(t,x))a;, 

^[A.(3)] ^ {D + X)q^^'^^^\t,x,x,x)-{D + X){£,{t,x))x. (2.18) 

Expanding the invariance condition, we find 

+ q(t>x + = 0, (2.19) 
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where the infinitesimal prolongations are as given above. If we put A = 0, we can get 
the Lie prolongation formula. Solving the invariance condition (12.191) we can obtain 
the explicit forms of 77 and A. If the given ODE admits Lie point symmetries, then 
the A-symmetries can be derived without solving the A-prolongation condition. In 
this case the A-symmetries can be deduced from the relation 


^ DU 
Q 


( 2 . 20 ) 


where D is the total differential operator and Q = rj — ^x. The associated vector 
field is given hj V = ^. The method of finding the integrals from A-symmetries 
can also be extended to the third-order ODEs as in the case of second-order ODEs. 


(/) Extended Prelle-Singer method 

In a series of papers Chandrasekar, Senthilvelan and Lakshmanan have devel¬ 
oped a stand-alone method, namely extended Prelle-Singer procedure, to investigate 
the integrability of the given ODE which may be of any order, including coupled 
ones. In the following, we recall briefly the extended Prelle-Singer procedure which 
is applicable for third-order ODEs (Chandrasekar et al. 2006). 

We assume that the ODE admits a first integral I{t,x,i,x) = C, with C 
constant on the solutions, so that the total differential gives 


dl = Itdt -I- Ij;dx + I±dx + Ixdx = 0 , ( 2 . 21 ) 

where the subscript denotes partial differentiation with respect to that variable. 
Rewriting equation in the form ^dt — dx = 0 and adding the null terms 

U{t,x,x,x)xdt — U(t,x,x,x)dx and S{t,x,x,x)xdt — S{t,x,x,x)dx to the latter, 
we obtain that on the solutions the 1 -form 

dt — Sdx — Udx — dx = 0. (2.22) 

Hence, on the solutions, the 1-forms (12.211) and (12.221) must be proportional, pro¬ 
vided (12.221) is a total differential. To ensure this we multiply (12.221) by the function 
R{t, X, X, x) which acts as the integrating factor for (12.221) so that we have on the 
solutions that 


— + Sx -\-Ux 


dl = R{cj) + Sx + Ux)dt — RSdx — RUdx — Rdx = 0. (2.23) 

Comparing now equations (j2.2ip with (12.231) we end up with a set of four relations 
which relates the integral(/), integrating factor(i?) and the null term(S'): 

= R{(j) -p iS' -p Ux), Ix = -RS, Ix = -RU, Ix = -R. (2.24) 

In order to determine S, U and R, we impose the compatibility conditions Itx = 
d-xt^ — ^xt: Rx — Ixt-, d^xx — dxx^ dxx — dxx and Ixx — dxx- Then we obtain the 
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following determining equations, 


D[S] = 

—(j>x + S(j)x + US, 

(2.25) 

D[U] = 

—(f^x U(j)x — s , 

(2.26) 

D[R] = 

-R{U + ^x), 

(2.27) 

Rx = 

RxS + RSx, 

(2.28) 

RxS = 

-SxR + RxU + RUx, 

(2.29) 

Rx = 

RxU -\- RUx, 

(2.30) 


where D = ^ + Solving the above system of over determined 

equations we can obtain the unknown functions S', U and R. From the known 
expressions, S, U and R, we can determine the integrals which appear on the left 
hand side of Ea. (l2.24l) by straightforward integration. 

Integrating Eg. (12.2411 we find 




I{t,x,x,x) = ri-r 2 -y 

■ / H + s'"' - "'ll''"' 


dx,(2.31) 


where 


ri 


R{(j) + Sx + Ux)dt, r 2 




We may note that every independent set (S, U, R) in (12.3111 defines a first inte¬ 
gral. 


3. Interconnections 


In the previous section we have discussed six specific analytical methods which are 
used to derive integrating factors, symmetries of various kinds, null forms and inte¬ 
grals associated with the third-order ODEs. A question which we raise now is what 
is the interconnection between these various methods, that is given any one of the 
quantities in the family, say multipliers, Darboux polynomials. Lie point symme¬ 
tries, adjoint-symmetries, A-symmetries or integrating factors and null forms, can 
one obtain the rest of the quantities in this family. To answer this question we can 
explore the hidden interconnections that exist between these functions and inter¬ 
link all these methods. To achieve this task we introduce certain transformations in 
the extended Prelle-Singer procedure which in turn connect globally all the above 
mentioned quantities. The details are given below (It may be noted that one can 
establish the same interconnection by taking any one of the other methods as the 
starting point). 

With the aid of the following transformations on the null forms of the PS 
method, see Ea. (l2.25L 


U=- 


D[V] 

V 


and 



(3.1) 


where V(t,x,x,x) and X{t,x,x,x) are two unknown functions and D is the total 
differential operator, Eos. (12.261) and (12.2511 respectively can be rewritten in the form 


D^[V] = D[V]<P^+V<Pi+X, 


(3.2) 
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and 

D[X] = (3.3) 

We introduce another transformation on the integrating factor as 

R = V/F, (3.4) 

where F{t,x,x,x) is a function to be determined, in (12.271) so that the latter can 
be rewritten in a compact form in the new variable F as 

D[F] = 4>iF. (3.5) 

One may note that Ea. (l3.5p is nothing but the Darboux polynomial determining 
equation ( 12 . 101 ) with the cofactor g = cj)^. 

We mention here that once we know the functions U and S the integrating 
factor R can be determined within the Prelle-Singer procedure itself. But to connect 
the integrating factors to Darboux polynomials these transformations are essential. 
More importantly the connection between the null forms and A-symmetries can be 
unearthed through the function V which appears in the transformations (EB and 
(133), as we see below. 

(a) Connection between symmetries and null forms 

Now we investigate how these expressions, namely U and S are interconnected 
with A-symmetries. In the case of second-order ODEs the A-symmetry is nothing 
but the null form with a negative sign (Muriel & Romero, 2009). This one-to-one 
correspondence came from the result that the S'-determining equation in the PS 
procedure differs only by a negative sign from that of the A-symmetry determining 
equation (Mohanasubha et al. 2014). However, in the case of third-order ODEs we 
have two null forms S and U which have to be connected to a single function A as 
demonstrated below. 

Let I{t, X, X, x) be a first integral of (12.11) then R = —I^ is an integrating factor, 
see Eg. (12.241) . The total derivative % = 0 gives 

R4> = It + xlx + xh- (3.6) 

Let I{t, X, X, x) also be a first integral of for some function A(t, x, x, x), then 

pD42)]j = o, (3.7) 

where pD’O)] jg given by -|- rj-^ + ^ -|- Here and 

are the first and second A- prolongations, respectively, whose explicit expressions 
are given in the first two expressions of (12.181) . Expanding Eq. (13.71) we get the 
following expression for A corresponding to V = ^-^ + rj-^, 

f)T 

+2X{'q^^^ - ^x) + X'^ig - ^x))—= 0, (3.8) 

where and 77 *-^) are the first and second Lie point prolongations which are given 
in Eg. (12.131) . For the vector field V = -^ the above expression (13.81) reduces to 

— (D[A] + y?)Ix = Ix F Xlx- (3.9) 
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Now we connect this expression which comes out from the A-symmetry analysis 
with the null forms in the extended PS procedure. In this regard we can deduce the 
following expressions from the last three equations of (I2.24L 

U=^, h^o. (3.10) 

Using the above in Eq. (13.91) we can obtain an equation which connects the null 
forms U and S with A in the form 


U = - 


(D[A] +X^ + S) 
A 


(3.11) 


Eq. (I3.11|) is essentially a redefinition of the A-function determining equation (13.911 . 
Thus in the case of third-order ODEs the null forms, U and S, are connected with 
the A-symmetries through a differential relation. This interconnection is brought 
out for the first time in the literature. We mention here that while deriving the 
relation (13.111) we assumed that X ^ 0. When A = 0 we have Ix = 0 (vide Eg. (13.Dll ') 
and in this case one of the null forms (S) vanishes. This result is also consistent 
with the extended PS procedure. 

For practical purpose we can rewrite the relation (13.1111 in terms of a single 
variable, say in V, as follows. Using dSS we can express X in terms of V and 
substituting the latter in the second expression in (EB we find S in terms of V. 
The resultant expression reads 


S=^{D^[V]-cl)xD[V]-cl)xV). 


(3.12) 


Now replacing the null functions, S and U, which appear in (13.1111 by V, we find 


D^[V] - {<Px + X)D[V] + (U[A] + A^ - <l,x)V = 0. (3.13) 


Substituting the known expressions (f>x and cj)x and A in (13.1311 and solving the 
resultant equation we can obtain V which in turn unambiguously fixes the null 
forms U and S through the relations given in Eg. (13.II) . In other words one can get 
the null forms U and S from A by finding the function V also. In this sense Eo. (l3.13p 
may also treated as a “second bridge” which connects A-symmetries with null forms. 
If we already know the null forms U and S, Eo. (l3.1ip yields the A-symmetries in a 
straightforward manner. In this sense one can determine (i) A from U and S and (ii) 
U and S from A. The expressions (13.111) and (13.131) interconnect the PS procedure 
and the A symmetry analysis. 


(6) Connection between Lie point symmetries and null forms 

The relation between Lie point symmetries and A-symmetries has already been 
established earlier (vide Eg. (12. 201) 1. Substituting this in (13.lip we can obtain an 
expression that relates Lie point symmetries with the null forms in the following 
manner: 

D^[Q] + UD[Q]+SQ = 0, (3.14) 

where Q = p — if is the characteristics. 

For the sake of completeness, in the following, we recall the following intercon¬ 
nection that is already known in the literature. 
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(c) Connection between Jacobi last multiplier and Darboux 
polynomials/integrating factors 

Using the expression (El we can relate the Darboux polynomials with the 
integrating factor and in fact the denominator of (Isa is nothing but the Darboux 
polynomials. 

By comparing Eas. (|2.10l) and (12.91) we can relate the Darboux polynomials with 
Jacobi last multiplier as 

M = F-\ (3.15) 

Using this relation, we can find the Jacobi last multiplier from the Darboux polyno¬ 
mials. So the integrating factor R in the PS method is the product of the function 
V and the Jacobi last multiplier M, that is R = VM. 

(d) Connection between adjoint-symmetries and integrating factors 

We rewrite Eas. (l2.25L (12.261) and (12.271) as a single equation in one variable, 
for example in R. Then the resultant equation reads 

D^[R] + D^[(f^R] - D[(j)iR] + (j),R = 0. (3.16) 

Comparing the adjoint of the linearized symmetry condition equation (|2.15|) and 
(13.161) one can conclude that the integrating factor R is nothing but the adjoint- 
symmetry A, that is 

R = A. (3.17) 

Thus the integrating factor turns out to be the adjoint-symmetry of the given 
third-order nonlinear ODE. 


(e) Connection between Lie point symmetries and Jacobi last multipliers 


The connection between Lie point symmetries and JLM is known for a long time 
in the form 


where 


A = 


1 

a 

6 


X 

X 

X 


(1) 

(2) 

Vi 

Wi 

Vi 




m 

V2 

V2 


(1) 

(2) 


m 

V3 


(3.18) 


(3.19) 


where (^i, ??i), {^ 2 ,^ 2 ) and {^ 3 ,^ 3 ) are three sets of Lie point symmetries (see below) 

(i) (i) (i) 

of the third-order ODE, ry) , ?72 and , i = 1, 2, are their corresponding first and 
second prolongations, respectively, and the inverse of A becomes the multiplier of 
the given equation (Nucci, 2005). 


(/) Comparison between interconnections for second- and third-order ODEs 

In the above said interconnections, some of the interconnections are common 
to both second- and third-order nonlinear ODEs except for their orders, while the 
others are different. Such common and differing connections are listed below. 
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1) Common features: 

The common connections are the ones between (i) Lie point symmetries and A- 
symmetries, (ii) Lie point symmetries and Jacobi last multiplier, (iii) adjoint- 
symmetries and integrating factors, and (iv) Darboux polynomials and Jacobi last 
multiplier. 

2) Differing features: 

(i) The uncommon relation is the connection between A-symmetries and null forms. 
In the cases of second- and third-order ODEs the connection between A-symmetries 
and null forms are entirely different. In the case of second-order ODEs the A- 
symmetry is nothing but the null form with a negative sign. This one-to-one cor¬ 
respondence came from the result that the S'-determining equation in the PS pro¬ 
cedure differs only by a negative sign from that of the A-symmetry determining 
equation (Muriel & Romero, 2009). In other words there is an one-to-one corre¬ 
spondence between the A-symmetries and null form S. However, in the case of 
third-order ODEs, we have two null forms {S and U) which have to be connected 
to the single function A. Our analysis shows that these two null forms are connected 
with the A-symmetry through a single expression. 

(ii) The connection between characteristics and null forms is also different in the 
case of second- and third-order ODEs. In the case of second-order ODEs, there 
exists only one null form which is directly connected with the characteristics, while 
in the case of third-order ODEs, the equation which connects the characteristics 
and null forms contains both the null forms. 


{g) Flow chart of the interconnections 

The above interconnections are clearly depicted in Fig.I. This may be compared 
with that of the second-order ODEs given as Fig.I in Mohanasubha et al. (2014). 


4. Illustration 


(a) Example 1: 

We consider the following third-order nonlinear ODE, namely 


XX 

X - = 0. 

X X 


(4.1) 


Equation grj is a sub-case of the general form of a scalar third-order ODE which 
is invariant under the generators of time translation and rescaling (Feix et al. 1997; 
Polyanin & Zaitsev 2003). A sub-case of equation (14.11) has been considered by both 
Bocharov et al. (1993) and Ibragimov & Meleshko (2005) in the form 'x — = 

0. They showed that it can be linearized to a linear third-order ODE through a 
contact transformation. The equation 'x — ^ = 0 has been considered by Euler 
& Euler (2004), who showed that it can be linearized to a free particle equation 
through a Sundman transformation. Equation (HU) was studied by Chandrasekar 
et al. (2005) from the integrability point of view using Prelle-Singer procedure. 
Here we consider Eq. (HU) from the perspective of deriving multipliers, Darboux 
polynomials. Lie point symmetries, adjoint-symmetries, A-symmetries, null forms. 


Article submitted to Royal Society 


13 


Lie point symmetries 



Jacobi last 
multiplier 


Darboux polynomials 


Adjoint symmetries 


Figure 1. Flow chart connecting Prelle-Singer procedure with other methods for 

third-order ODEs 


integrating factors, sequentially and demonstrate the effectiveness of exploring the 
interconnections. From these quantities we also derive integrals and the general 
solution of this equation for the sake of completeness. 

We begin our analysis with Lie point symmetries. Equation 63) admits a set 
of three parameter Lie point symmetries of the form 




at’ 


0,2 = —3 


dx' 


— t 


df 


(4.2) 


The characteristics associated with the above vector fields is given by 


Qi = -i, Q 2 = -X, Qs = -tx. 


(4.3) 


Using the relation A = —^ we can derive the A-symmetries associated with the 
above vector fields which in turn read 


Ai 




X + tx 
tx 


(4.4) 


Once we know the A-symmetry we can find the null forms U and S. To obtain them 
we determine the function V using the relation (I3.13p . We obtain 


Vi = X, V 2 = x^ — XX, V 3 = tx^ — x(x + tx). (4.5) 


We can fix the allowed forms of the function X from the known expression of V 
through the relation (|3.2I) . The resultant expressions read 


Xi = 


XX 


X2=- 


• 2 •• 

XX 


X3 = - 


xx(x -|- tcc) 


(4.6) 
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Using the functions X and V we can obtain the null forms U and S respectively, 
vide Eq. (EB , as 


) ^‘2 / • 2 " \ ’ 
X x(x^ — XX) 


S3 = 


xx(x + tcc) 


and 


C/r=-^, 

X 


U 2 = 


XX 


U3 = - 


x{—tx^ + x(x + tx))' 
xx( 2 x + tx) 


(4.7) 


(4.8) 


x(x^ — xx)^ x(tx^ — x{x + tx)) 

Substituting the null forms U and S in Ea. (l2.27p and solving the resultant equation 
we find the integrating factors for dUD in the form 




-1 


2 x 


7?2 .2 5 

XX X‘^ X 


i?3 = 


tx^ — x(x + tx) 


2 xx 


x{2x'^ —xx) 


(4.9) 


The functions Ri , R 2 and R 3 also become adjoint-symmetries of as well which 
can be directly verified by substituting them in (12.151) . Using the expressions for R 
and V we can fix the Darboux polynomials for the given equation by recalling the 
relation F = ^. The Darboux polynomials turn out to be 

• 2 

Ell =Ei 2 =Ei 3 = 2a;i:^. (4.10) 

The polynomials Eii,Ei 2 and E 13 share the common cofactor f + ^- From the 
property that the ratio of Darboux polynomials which share the same cofactor is 
the first integral, we can find another Darboux polynomial which shares the cofactor 
as the above. Doing so, we have found two more Darboux polynomials satisfying 
the equation D[E] = (| -|- ^)E as 


F 2 = x{2x^ — xx), E 3 = XX. 


(4.11) 


We can find more number of Darboux polynomials using the property of Darboux 
polynomial that the combination of Darboux polynomials is also a Darboux poly¬ 
nomial. By multiplying an integral with the Darboux polynomial we can obtain 
more Darboux polynomials. The ratios of the Darboux polynomials and lead 
to the first integrals /i and I 2 (vide Eg. (14.121) '). respectively. 

The JLM associated with the given equation can be obtained from the Darboux 
polynomials (14.101) by recalling the relation F = M~^. 

Once we know the null forms and integrating factors we can construct the asso¬ 
ciated integrals of motion by evaluating the integrals given in (I2.3ip . Our analysis 
shows that 


XX 


X x^ 


h — l 2 = x[ -ry , h = - -ZT\ -( 2 x^ - xx) + tan 


2x \ X 


x( 2 x^ — xx) 

(4.12) 

It is a straightforward matter to verify that Ii , I 2 and I 3 are the three independent 
integrals for the given Eq. 611). From the knowledge of /i, I 2 and I 3 we get the 
general solution in the form 


x(t) 




(4.13) 


The above said results are tabulated in Table 1. 
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( 6 ) Example 2: 


We consider another interesting example (Bluman & Anco, 2002) 


'x 


6tx^ 6x^ 

—rj- H-• 

X 


(4.14) 


The adjont symmetries of this equation were first worked out by Bluman and Anco 
(2002). Subsequently the extended Prelle-Singer method was applied to this equa¬ 
tion by Chandrasekar et al. (2005) and they have derived the null forms, integrating 
factors and integrals through this method. The exact expressions are given in Table 
2. The procedure given in the previous example may be followed for this example as 
well to derive the quantities displayed in Table 2. Here also because of nonunique¬ 
ness of S, U with A, one may find two independent integrals I 2 and from the 
same A-function. The vector field Hi provides one independent integral Ii and the 
vector field H 2 gives the two independent integrals I 2 and The vector field H 3 
provides only a dependent integral. 

Using the integrals Ii,l 2 and the solution of Ea. (l4.14p can be written down 
implicitly as 


3 t(/ 2 (/i - 2x) + 9thf + / 2 ((/i + ‘ixf - mhf 

-{h + 2 a;)(/ 2 (/i -b 2x) + %th){{h + 2xY - Uth) = 0. (4.15) 


(c) Example 3: 


Finally, we consider the following example. 


'x 


x(x — 1 ) 
a; -b 1 


Eg. (14. 161) admits two Lie point symmetries which are given by 


Hi 


'We 


H2-t^ + 2(x + l)^. 


(4.16) 


(4.17) 


Since Ea. (l4.16p admits only two point symmetries, we start our analysis by ex¬ 
ploring the null forms and the integrating factors in the PS method with the help 
of Eas. (|2.25ll . (I2.26|) and (12.2711 . Integrating these three equations, we obtain the 
following particular solutions for the null forms Si and Ui, i = 1, 2, 3, namely 


Si 

53 


\ — x a:^(i-b 1 )-b 2 a;(x-b 1 )-b 2 i 

x + \' ^ x{x + \){x + 2) 

{x — 1 ) {t {x^{x — 1 ) — (a; -b l)i^) -b (a; -b l)a;x) 

(a: -b 1 ) {t{x -b l)x‘^ -b x‘^{t — tx) -b (a; -b l)x{x — 2 )) ’ 


r/ = n TT = ~^(^ + 1)^ Tj ^ _ 2(a:-bl)(l-:r):r _ 

^ ^ a:(x-b 2 ) ’ ^ t{x + l)x‘^ + x^{t — tx) + {x+ \)x{x 

The associated integrating factors are found to be 


2 )' 

(4.18) 


Ri 


-1 

X -b 1 ’ 


i?2 


x{x + 2 ) 
2(x-bl)’ 


i?3 


—t{x + l)iP + tx^{x — 1 ) — (x -b l)x(x — 2 ) 
2 (a; -b 1 )^^ ((a; -b l)x'^ - x^{x - 1 )) 

(4.19) 
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Once we know the null forms S and U, we can deduce the A-symmetries by solving 
Eg. (13.111) . Doing so, we find 


Ai — A2 



x{x — 1) 
{x + l)x' 


(4.20) 


The null form pair { 83 ^ U 3 ) provides another A-symmetry which is given by 


A 2 


X — tx 
2 + 2x — tx 


(4.21) 


The A-functions Ai(= A 2 ) and A 2 can also be obtained directly from the Lie point 
symmetries (14.171) . In this example also, we observe that for a single A one can have 
multiple null forms. From the null forms, the functions V and X can be determined 
with the help of Eq. (EB as 


El 


^1 

^3 


1 , H = x{x + 2), Vi = + + + !)»(» 

X +1 

1 — x x^ {x + 1 ) -\-2x{x-\-1)-\-2x 

{x — 1) {t {x^{x — 1) — (a; + l)i^) + (a: + l)xx) 

{x + iy 


2 ) 


(4.22) 


Since the functions V and R are known now, the Darboux polynomials admitted 
by Eg. (14.161) are found by exploiting the relation (13.41) . We observe 


Fi = -(ai + l), F 2 = 2 (a: + 1 ), F 3 =-2 


\ — X 

a; + 1 


((x + l)ai^ — a;^(ai — 1)) . (4.23) 


The JLM associated with the given equation can be obtained from the Darboux 
polynomials (I4.23|) by recalling the relation F = M~^. Jaboci last multipliers are 
given by 


Ml = - 


1 + X 


M 2 = 


1 


2 x 2 


M3 = -- 


((a; + l)i^ - x2(x - 1)) 


(4.24) 


Once we know the null forms and integrating factors we can construct the associated 
integrals of motion by evaluating the integrals given in (12.311) . Our analysis shows 
that 


h 

h 


X + X 

X + 1 
tan“^ 


h = 


x^ — x^(x + 1) + X (x^ — 2x) 


2 (x + l) 




i:(x — 1 ) 


- t 


1 — x 
X + 1 


(4.25) 


Using the above three integrals, we can derive the general solution of (|4.16l) in the 
following form 


c(i) = 


h - I sin(/3 + y/l - lit) 
"^1 ^ 


(4.26) 


where I = y/if + 2/2 — 2 / 1 / 2 . 


Article submitted to Royal Society 

































17 


5. Panorama of interconnections 

In the previous section, starting from Lie point symmetries we derived all other 
quantities. Since the interconnections are global one can consider any other quantity 
in this family, and derive the rest of them, see the connection diagram (Fig.l). In 
this section, we consider two such cases, namely (i) integrating factors as starters 
and (ii) Darboux polynomials as starters and demonstrate the method of deriving 
all other quantities from them. 

(a) From integrating factors to others 

In this subsection we demonstrate that starting from the integrating factor 
we can derive all other quantities for the example (ED. Suppose an integrating 
factor Ri (i?i = — ^) is given as a starter (vide Eg. (14.91) 1. Then we can get the 
corresponding null form Ui from Ea. (l2.27p which exactly matches with the first 
expression given in Eg. (14.81) . Substituting C/i in (El and rewriting the resultant 
equation we obtain 

7?[Pi]- =0. (5.1) 

A particular solution to this equation is Vi = — i which in turn agrees with (14. 5p . 
Using Vi we can get an expression for Xi through the relation Eq. (El which also 
exactly matches with the first expression given in (14.61) . Once we have Xi and Vi 
we can obtain the second null form Si with the help of Eq. (E- Once we know Si 
and Ui we can construct the A-symmetry through the expression (13.1111 . that is 

D[X]+X‘^ -^X-- =0. (5.2) 

X X 

A particular solution to this equation is A = j which also agrees with the first 
expression given in (lOl) . The corresponding A-symmetry is ^. Using the quantities 
i?i and Vi in E we can find the Darboux polynomials as given in the first 
expression in (I4.10|) . In this way we can derive the rest of the quantities from an 
integrating factor. The exact expressions for the other quantities can be found in 
Table 1. The associated integral turns out to be h as expected. The procedure is 
exactly the same for the other two integrating factors i ?2 and i ?3 which are given 
in (E. The integrating factors i ?2 and i ?3 provide exactly the same expressions 
{S 2 ,U 2 , R 2 ) and {Ss,!!^, R 3 ). We also note here that if we start from some other 
integrating factor other than the above three then proceeding in the same manner 
as outlined above one can get their associated null forms and Darboux polynomials. 
However these expressions do not lead to any new integrals. 

(&) From Darboux polynomials to others 

Suppose Darboux polynomials which share the same cofactor are given as first 
information and we have to determine the rest of the quantities. Erom the Darboux 
polynomials we can find the integrals using the ratios. We assume that the two 
Darboux polynomials F 12 = xx^ and F 3 = xx, vide Eos. (14.101) and (14.111) . are 
given. Then the ratio gives the integral h = (vide Ea. ()4.12[) '). Once we know 
the integral we can obtain the integrating factor readily by evaluating the last 
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expression given in (I2.24p which exactly matches with the form of Ri given in 
(gH). The method of deriving the rest of the quantities from i?i is outlined in the 
previous sub-section. In this way, from Darboux polynomials one can derive the rest 
of the quantities. The procedure is the same for the other Darboux polynomials. 

One may assume that the last multipliers are given and then consider the task 
of determining the other quantities. Using the relation F = M~^ we can find the 
Darboux polynomials. Once Darboux polynomials are known the procedure given 
in Sec.4(c) may be followed to derive the rest of the quantities. 


6. Conclusion 

In this work, we have demonstrated the interconnection between Prelle-Singer 
method (or any one of the methods studied in this paper as the starting point) 
with the other existing well known methods in the literature such as Jacobi last 
multipliers, Darboux polynomials. Lie point symmetries, adjoint-symmetries and A- 
symmetries in the case of third-order nonlinear ODEs. For this purpose we started 
with the PS method. In the PS method, the quantities, namely (i) null forms U and 
S and (ii) integrating factors R which are determined by three first-order equations 
play a major role. By introducing suitable transformations to the null forms U and 
S and to the integrating factor R, we have related these three quantities with the 
above other quantities. While relating the PS method with Darboux polynomials 
and Jacobi last multiplier, we introduced the transformation ^ in the integrat¬ 
ing factor equation. To relate the adjoint-symmetries with the integrating factor 
we rewrite the three first-order equations in the PS method in terms of a single 
variable R. We then demonstrate that this third-order equation in the variable 
R coincides with the adjoint-symmetry equation. The difficult and unknown con¬ 
nection between A-symmetries and the null forms has been brought out by using 
the compatibility between the A-determining equation and the determining equa¬ 
tion for the null forms in the PS method. We have recalled the known connection 
A = to relate the Lie point symmetries and the null forms. In this way we 
have established the interconnections between the null forms, integrating factors, 
adjoint-symmetries, A-symmetries, Lie point symmetries, Darboux polynomials and 
Jacobi last multipliers. We have observed that some of the connections are common 
for both second-order as well as third-order ODEs, while the others are specifically 
applicable to third-order ODEs. We have illustrated these interconnections with 
three definitive examples discussed in the literature. Currently we are extending 
the above said procedure to n*^-order nonlinear ODEs. We have obtained some 
interesting results in this direction. The results will be published in the near future. 
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A 

V 

u 

X 

s 

F 

R 

I 

dt 

X 

X 

X 

X 

X 

XX 

X 

X 

X 

—xx^ 

-1 

XX 

XX. 



x^x 

1 2x 1 X j 

\ X ' X j 

XXX 

X 

X 

xx^ 

X 

X 

II 

ox 

X 

X^ — XX 

XX^ 

x'^x 

x^x 

— xx^ 

2x 2 

= 

X 

x(x^ —xx) 

X 

x{x^ — xx) 

2 

X^ X 

^ dt 

x-\-tx 

tx^ — x(x + tx) 

xx(2x-i-tx) 

xx(x-\-tx) 

xx{x-\-tx) 

0 ™™, / 4{2±3_a:i;) 

tx^ —x{x-\-tx) 


tx 

x(tx'-^ —x(x-i-tx)) 

X 

x{ — tx^-\-x{x-\-tx)) 

" y X 









~XX Y ^ 

-|-tan“^ ,/ ,„.a— tt^x 

y x{2x'^ — xx) 


Table 1. Vector fields (n), A-symmetries, null forms {U,S), Darboux polynomials {F), integrating factors (R) and integrals (!) admitted by 
Eq. (HU, along with the quantities V and X defined by Eqs. dUl) and (iTSll . respectively. 


n 


' dt 


x+tx 


V 


u 


2xx+6tx'^ 


3xx-\-6tx^ 


Axx-\-&tx‘‘ 


X 


( x-^2tx 


R 


h = 6tx — 2x + 


3^3 2t+f 


I2 = 3tx^ + ^ 


h = V 




Table 2. Vector fields (17), A-symmetries, null forms {U,S), Darboux polynomials {F), integrating factors (R) and integrals (/) admitted by 
Eq. (14.1411 . along with the quantities V and X defined by Eos. (13.211 and (13.311 . respectively 


o 

(M 
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